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Abstract. We consider a coupled bistable A^-particle system on R N driven by a 
Brownian noise, with a strong coupling corresponding to the synchronised regime. 
Our aim is to obtain sharp estimates on the metastable transition times between the 
two stable states, both for fixed N and in the limit when N tends to infinity, with 
error estimates uniform in N. These estimates are a main step towards a rigorous 
understanding of the metastable behavior of infinite dimensional systems, such as 
the stochastically perturbed Ginzburg-Landau equation. Our results are based on 
the potential theoretic approach to metastability. 
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1. INTRODUCTION 

The aim of this paper is to analyze the behavior of metastable transition times for 
a gradient diffusion model, independently of the dimension. Our method is based 
on potential theory and requires the existence of a reversible invariant probability 
measure. This measure exists for Brownian driven diffusions with gradient drift. 

To be specific, we consider here a model of a chain of coupled particles in a 
double well potential driven by Brownian noise (see e.g. [2]). I.e., we consider 
the system of stochastic differential equations 

dX e (t) = -VF^ N (X e (t))dt + V2~edB(t) } (1.1) 

where X e (t) e R N and 

ieA ^ ' ieA 

with A = Z/iVZ and 7 > is a parameter. B is a N dimensional Brownian motion 
and e > is the intensity of the noise. Each component (particle) of this system is 
subject to force derived from a bistable potential. The components of the system 
are coupled to their nearest neighbor with intensity 7 and perturbed by indepen- 
dent noises of constant variance e. While the system without noise, i.e. e = 0, has 
several stable fixpoints, for e > transitions between these fixpoints will occur at 
suitable timescales. Such a situation is called metastability. 

For fixed N and small e, this problem has been widely studied in the literature 
and we refer to the books by Freidlin and Wentzell Q9D and Olivieri and Vares [TT5T1 
for further discussions. In recent years, the potential theoretic approach, initiated 
by Bovier, Eckhoff, Gayrard, and Klein |5]] (see O for a review), has allowed to 
give very precise results on such transition times and notably led to a proof of 
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the so-called Eyring-Kramers formula which provides sharp asymptotics for these 
transition times, for any fixed dimension. However, the results obtained in H5U do 
not include control of the error terms that are uniform in the dimension of the 
system. 

Our aim in this paper is to obtain such uniform estimates. These estimates con- 
stitute a the main step towards a rigorous understanding of the metastable be- 
havior of infinite dimensional systems, i.e. stochastic partial differential equations 
(SPDE) such as the stochastically perturbed Ginzburg-Landau equation. Indeed, 
the deterministic part of the system (11.1D can be seen as the discretization of the 
drift part of this SPDE, as has been noticed e.g. in [3]. For a heuristic discussion 
of the metastable behavior of this SPDE, see e.g. lfT3H and IfTTH - Rigorous results 
on the level of the large deviation asymptotics were obtained e.g. by Faris and 
Jona-Lasinio iTTOll . Martinelli et al. |fT4lh and Brassesco [0]. 

In the present paper we consider only the simplest situation, the so-called syn- 
chronization regime, where the coupling 7 between the particles is so strong that 
there are only three relevant critical points of the potential F ljN (11. 2D . A gen- 
eralization to more complex situations is however possible and will be treated 
elsewhere. 

The remainder of this paper is organized as follows. In Section 2 we recall briefly 
the main results from the potential theoretic approach, we recall the key properties 
of the potential F 7 at, and we state the results on metastability that follow from the 
results of fl5j for fixed N. In Section 3 we deal with the case when iV tends to 
infinity and state our main result, Theorem 13.11 In Section 4 we prove the main 
theorem through sharp estimates on the relevant capacities. 

In the remainder of the paper we adopt the following notations: 

• for t & K., L^J denotes the unique integer k such that k < t < k + 1; 

• t d = inf{t > : X t e D] is the hitting time of the set D for the process 

(**); 

• B r (x) is the ball of radius r > and center x e W N ; 

• ior p > 1, and (xk)k=i a sequence, we denote the L p -norm of x by 

/ N \ VP 
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2. Preliminaries 

2.1. Key formulas from the potential theory approach. We recall briefly the 
basic formulas from potential theory that we will need here. The diffusion X e is 
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the one introduced in (11. ID and its infinitesimal generator is denoted by L. Note 
that L is the closure of the operator 

L = ee F ^ N/£ Ve~ F ^ N/£ V. (2.1) 

For A, D regular open subsets of M. N , let h A>D (x) be the harmonic function (with 
respect to the generator L) with boundary conditions 1 in A and in D. Then, for 
x G (A U D) c , one has h AtD (x) = ^x[t~a < Td]- The equilibrium measure, e A> o, is 
then defined (see e.g. fl8D) as the unique measure on dA such that 

h A , D {x)= f e- F ^/ £ G D o(x,y)e A , D {dy), (2.2) 

JdA 

where G D c is the Green function associated with the generator L on the domain 
D c . This yields readily the following formula for the hitting time of D (see. e.g. 
0): 

/ E z [r D ]e- F ^/ £ e A!D (dz) = [ h A , D {y)e- F ^' £ dy. (2.3) 

JdA JD C 

The capacity, cap(A, D), is defined as 

cap(A,D)= [ e- F ^ N(z)/£ e AD (dz). (2.4) 

JdA 

Therefore, 



e 



Mz)/£ e A , D (dz) 



cap(A, D) 

is a probability measure on dA, that we may call the equilibrium probability. The 
equation (12. 3D then reads 

E z [T D ]z/ A)D (ck) = E^ D [r D ] = — — . (2.6) 

9j4 cap^/i,^; 

The strength of this formula comes from the fact that the capacity has an alternative 
representation through the Dirichlet variational principle (see e.g. UTTJ]), 

cap(A,D) = inf (2.7) 

where 

U= \h G W 1 ' 2 ^ ,e- F ^ N{u)/£ du) \Vz,h{z) G [0, 1] , h ]A = 1 , h\ D = o}, (2.8) 
and the Dirichlet form $ is given, for h eH, as 

$(h)=e f e- F ^ N{u)/£ \\Vh(u)\\ 2 2 du. (2.9) 

J(AUD) C 

Remark. Formula (I2.6D gives an average of the mean transition time with respect 
to the equilibrium measure, that we will extensively use in what follows. A way 
to obtain the quantity E 2 [t d ] consists in using Holder and Harnack estimates Ifl2ll 
(as developed in Corollary I2.3D [51. but it is far from obvious whether this can be 
extended to give estimates that are uniform in N. 

Formula (I2.6D highlights the two terms for which we will prove uniform esti- 
mates: the capacity fProposition l4.3D and the mass of h A> o (Proposition 14. 9D . 
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2.2. Description of the Potential. Let us describe in detail the potential F ltN , its 
stationary points, and in particular the minima and the 1 -saddle points, through 
which the transitions occur. 

The coupling strength 7 specifies the geometry of F 7jA r. For instance, if we set 
7 = 0, we get a set of N bistable independent particles, thus the stationary points 
are 

£* = (&,..., Cat) Vie[l,iV], & e {-1,0,1}. (2.10) 
To characterize their stability, we have to look to their Hessian matrix whose signs 
of the eigenvalues give us the index saddle of the point. It can be easily shown 
that, for 7 = 0, the minima are those of the form (I2.10D with no zero coordinates 
and the 1-saddle points have just one zero coordinate. As 7 increases, the structure 
of the potential evolves and the number of stationary points decreases from 3^ to 
3. We notice that, for all 7, the points 

i± = ±(l,l,...,l) O = (0,0,...,0) (2.11) 

are stationary, furthermore I± are minima. If we calculate the Hessian at the point 
O, we have 



/-I +7 - 



V 2 R 



7,JV 



(O) 



V 



1 
2 

1 + 7 

2 















\ 







(2.12) 



-1 -1 + 7/ 



whose eigenvalues are, for all 7 > and for < k < N — 1, 



A 



k,N 



(2.13) 



Set, for k >1,% 



N 



2 sin 2 (kn/N) 
Afc,AT = \N-k,N = 

-Vat = Aq = — 1. 



Then these eigenvalues can be written in the form 

1 



-1 + 4-, 1< k < N 



(2.14) 



Note that (7^)^=1 is a decreasing sequence, and so as 7 increases, the number 



of non-positive eigenvalues (Xk,N) k=0 decreases. When 7 > 7(7 the only negative 
eigenvalue is — 1 . Thus 

1 



N 



7i 



2sin 2 (7r/iV) 



(2.15) 



is the threshold of the synchronization regime. 



Lemma 2.1 (Synchronization Regime). If 7 > 7^7 the only stationary points of F J:N 
are I± and O. I± are minima, O is a 1 -saddle. 

This lemma was proven in H2D by using a Lyapunov function. This configuration 
is called the synchronization regime because the coupling between the particles is 
so strong that they all pass simultaneously through their respective saddle points 
in a transition between the stable equilibria (I±). 

In this paper, we will focus on this regime. 
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2.3. Results for fixed N. Let p > and set B± = B P (I±), where B p (x) denotes 
the ball of radius p centered at x. Equation (12.61) gives, with A = B_ and D = B + , 

f B c h B _, B+ (y)e- F -<My)/zdy 

E — + W ^ + cap(B-,B + ) ' (2 - 16) 

First, we obtain a sharp estimate for this transition time for fixed N: 

Theorem 2.2. Let N > 2 be given. For 7 > 7 f = — Zet v 7 ^ > p > e > 0. 

Then 

E, B _, fl+ [r B+ ] = 27rc w e^(l + 0(v/£|ln£|3)) (2.17) 

with 

c^v= TT 

L 1 'J fc=1 
where e(N) = 1 if N is even and if N is odd. 

Remark. The power 3 at In e is missing in HSJ by mistake. 

Remark. As mentioned above, for any fixed dimension, we can replace the prob- 
ability measure v b _,b + by the Dirac measure on the single point / , using Holder 
and Harnack inequalities Q5J . This gives the following corollary: 

Corollary 2.3. Under the assumptions of Theorem \2.2\ there exists a > such that 



9 + -X 



(2.18) 



E/_[r B+ ] = 2irc N eK(l + 0(y/e\ lne| 3 ). (2.19) 

Proof of the theorem. We apply Theorem 3.2 in 0. For 7 > 7^ = 2 S m a Twjv j j ^ et us 
recall that there are only three stationary points: two minima I± and one saddle 
point O. One easily checks that F JjN satisfies the following assumptions: 

• F 7i 7v is polynomial in the (xj) ie A and so clearly C 3 on WL N . 

x— >oo 

• ||VF 7>A r(x)|| 2 ~ as ||x|| 2 ->■ 00. 

• As Ai 7iW (i) ~ 3||x||| (||x|| 2 -> 00), then || VF 7>Af (x)|| - 2AF 7)JV (i) ~ ||x||g. 
The Hessian matrix at the minima I± has the form 

V 2 F 7i7V (/±) = V 2 F 7i7V (0) + 3Id, (2.20) 
whose eigenvalues are simply 

Vk,N = Afc jv + 3. (2.21) 



Then Theorem 3.1 of [5] can be applied and yields, for y/N > p > e > 0, (recall 
the the negative eigenvalue of the Hessian at O is —1) 



*b_,b + [ t b+] = — 7 =__== — (1 + (v ^|l Ile | )). (2.22) 



Finally (I2TT4D and (EOTD give: 

I JV-l I I JV-1 I 

det(V 2 F 7 ^(/_))=n^JV = 24^ J] ^ = 2^(1 + 7)^ U l + ^v) 



fc=0 fe=l fe=l 



(2.23) 
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N-1 



|det(V 2 F 7iiV (0))|= n^ = A^ f[ A 2 . ^ = (2 T - l) e W f[ (l 

k=0 k=l k=l ^ 

Then, 



_7_ 



c N 



v /det(V 2 F 7 , jV (/_))" 
7|det(V 2 J P 7)JV (0))| 

and Theorem 12.21 is proved. 



2 + 2 7 



n 

k=l 



2 4- -2- 



(2.24) 

(2.25) 

□ 



Let us point out that the use of these estimates is a major obstacle to obtain a 
mean transition time starting from a single stable point with uniform error terms. 
That is the reason why we have introduced the equilibrium probability. However, 
there are still several difficulties to be overcome if we want to pass to the limit 

iV t oo. 

(i) We must show that the pref actor c N has a limit as N | oo. 

(ii) The exponential term in the hitting time tends to infinity with N. This 
suggests that one needs to rescale the potential F 7 jv by a factor 1/N, or 
equivalently to increase the noise strength by a factor N. 

(iii) One will need uniform control of error estimates in iV to be able to infer 
the metastable behavior of the infinite dimensional system. This will be the 
most subtle of the problems involved. 



3. Large N limit 

As mentioned above, in order to prove a limiting result as iV tends to infinity, we 
need to rescale the potential to eliminate the iV-dependence in the exponential. 
Thus henceforth we replace F ltN (x) by 

G y>N (x) = N~ 1 F j ^n(x). (3.1) 

This choice actually has a very nice side effect. Namely, as we always want to 
be in the regime where 7 ~ 7^ ~ iV 2 , it is natural to parametrize the coupling 
constant with a fixed /j, > 1 as 

Then, if we replace the lattice by a lattice of spacing 1/N, i.e. (xi) ie A is the dis- 
cretization of a real function x on [0, 1] (x{ = x(i/N)), the resulting potential con- 
verges formally to 

G^(x) ^ jT (j[x{8)]' - i[*( S )] 2 ) ds + £- 2 £ [ ^^ds, (3.3) 
with x(0) = x(l). 

In the Euclidean norm, we have ||/±||2 — VN, which suggests to rescale the size 
of neighborhoods. We consider, for p > 0, the neighborhoods B± = B py r^(I±). The 
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volume V(B^) = V(B+) goes to if and only if p < \j2ixe, so given such a p, the 
balls B± are not as large as one might think. Let us also observe that 



i t 1 . .... 



x\\ 2 — > ||x||l2[o,i] = / \x(s)\ ds. (3.4) 







Therefore, if x e B+ for all N, we get in the limit, \\x — 1\\l 2 [o,i] < P- 

The main result of this paper is the following uniform version of Theorem 12.21 
with a rescaled potential G 7i jv- 

Theorem 3.1. Let p e]l, oo[, then there exists a constant, A, such that for all N > 2 
and all e > 0, 

-^E„ [t b jv] = 27r Cjv e 1/4£ (l + R(e,N)), (3.5) 



where c N is defined in Theorem \2.2\ and \R(e,N)\ < AWe\ lne\ 3 . In particular, 



where 



lim lim ^e~^E UBN gN [t b n] = 2nV(p) (3.6) 
^) = n^<«- C3.73 

fc=l r 

Remark. The appearance of the factor l/N may at first glance seem disturbing. It 
corresponds however to the appropriate time rescaling when scaling the spatial 
coordinates i to i/N in order to recover the pde limit. 

The proof of this theorem will be decomposed in two parts: 



• convergence of the sequence c N (Proposition [3T2]); 

• uniform control of the denominator (Proposition I4.3D and the numerator 
(Proposition g2D of Formula J235b . 

Convergence of the prefactor c N . Our first step will be to control the behavior of 
cat as iV t oo. We prove the following: 

Proposition 3.2. The sequence c N converges: for p > 1, we set 7 = p^i, then 

lim c N = V(p), (3.8) 

iVtoo 



with V(p) defined in (13771) . 

Remark. This proposition immediately leads to 

Corollary 3.3. For p e]l, oof, we set 7 = p^ , then 

_j_ 

^ m H ^ C ~W^ b^v I t b-} = ^V(p). (3.9) 

JVtoo e^O iV B - - B + + 

Of course such a result is unsatisfactory, since it does not tell us anything about 
a large system with specified fixed noise strength. To be able to interchange the 
limits regarding e and N, we need a uniform control on the error terms. 
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Proof of the proposition. The rescaling of the potential introduces a factor j* for the 
eigenvalues, so that (12.22D becomes 



E 



Ar-^/Vdet(V 2 F 7 ^(/_)) 
2irNc N e^(l + ©(v^llnel 3 )). 



e\ 3 )) 



(3.10) 



Then, with u N 



k — ,f > 



Cat 



1 - 



2 + 2/ryf 



N 



n 



fe=l 

N\N-1 



(3.11) 



To prove the convergence, let us consider the (7^ )^ =1 • For all fc > 1, we have 



21 

7 fc W sin 2 (£) 
Hence, u% „—>_ _v k = ^p. 



JV-,.+oo 



Thus, we want to show that 

+00 



1 

iV2 



(3.12) 



c N 



l[(l-v h )=VQi). 



Using that, for < t < f , 



fe=i 



f 2 



< t 2 (l - — ) < sin 2 (t) < t 2 , 



(3.13) 



(3.14) 



we get the following estimates for set a = ( 1 — ^ J, for 1 < < AT/2, 

k 2 



ak z 



k 2 7l 2 

3N 2 



< 



7T 
7? 



< 



Then, for N > 2 and for all 1 < k < N/2, 

k 4 Tr 2 >y N 

Let us introduce 



k 2 n 2 



A: 2 tt 2 



< 



(3.15) 



(3.16) 



m— 1 



n t^, ro = n (1 
fc=i /c=i 



Then 



In 



[7 



N 



< 



fc=i fe=i 

Using (137151 and (I3TT61) . we obtain, for all 1 < k < N/2, 



N 
k 



(3.17) 



(3.18) 



N 

v k -u k 


3/U 




-k 2 




1 - v k 


(-1 + /^ 2 ) 





4^2 



< 



Ilk 7T 



< 



C 



A^ 2 (-1 + /i£; 2 ) (2 + ^afc 2 ) ~ N 2 



(3.19) 
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with C a constant independent of k. Therefore, for N > N Q , 

C 



In - 




In ^ 


1 - v k 







1 + 



N 

v k - < 

1 - v k 



< 



N 2 ' 



Hence 



In 



U. 



N,N 



Vr 



N 



< 



N n- 



0. 



(3.20) 



(3.21) 



As J2\ v k\ < +°°! we g et nm Af->+oo Vn = V{n) > 0, and thus (I3.13D is proved. □ 



4. ESTIMATES ON CAPACITIES 

To prove Theorem 13.11 we prove uniform estimates of the denominator and 
numerator of (I2.6D . namely the capacity and the mass of the equilibrium potential. 

4.1. Uniform control in large dimensions for capacities. A crucial step is the 
control of the capacity. This will be done with the help of the Dirichlet principle 
A2.7D . We will obtain the asymptotics by using a Laplace-like method. The expo- 
nential factor in the integral (I2.9D is largely predominant at the points where h is 
likely to vary the most, that is around the saddle point O. Therefore we need some 
good estimates of the potential near O. 

4.1.1. Local Taylor approximation. This subsection is devoted to the quadratic ap- 
proximations of the potential which are quite subtle. We will make a change of 
basis in the neighborhood of the saddle point O that will diagonalize the quadratic 
part. 

Recall that the potential C7 7 >N is of the form 

G,A*) = [Id - D]x) + ^IMIt (4.1) 

where the operator D is given by D = 7 [id — ~(E + E*)l and (Hx)j = x j+ i. The 
linear operator (Id — D) = — V 2 F 1>N (0) has eigenvalues —Xk,N and eigenvectors 
Vk,N with components v^nU) = u jk , with u = e %2lT ^ N . 
Let us change coordinates by setting 

JV-l 

xj = tu~ jk x k . (4.2) 

Then the inverse transformation is given by 

1 JV " 1 

j=Q 

Note that the map x — > x maps R N to the set 

R N = {x e C N : x k = I^~ k } (4.4) 
endowed with the standard inner product on C N . 

Notice that, expressed in terms of the variables x, the potential (II. 2D takes the form 

1 N ^ 1 
G 7>N {x{x)) = — KN\x k \ 2 + ^\\x(x)\\t. (4.5) 

k=0 
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Our main concern will be the control of the non-quadratic term in the new coordi- 
nates. To that end, we introduce the following norms on Fourier space: 



N-l \ Vp 

I 1 



i V i 



N ^ 1 1 ) N 1 Ip 

i=0 / 

The factor 1/N is the suitable choice to make the map x — > x a bounded map 
between L p spaces. This implies that the following estimates hold (see IfToll . Vol. 
1, Theorem IX. 8): 

Lemma 4.1. With the norms defined above, we have 

(i) the Parseval identity, 

1Mb = \\xh,T, (4.7) 

and 

(ii) the Hausdorjf-Young inequalities: for 1 < q < 2 and p^ 1 + q^ 1 = 1, there 
exists a finite, N -independent constant C q such that 

\\x\\ p < C q \\x\\ qtF . (4.8) 

In particular 

\\x\\i < C4/ 3 ||s||4/ 3r 77. (4.9) 

Let us introduce the change of variables, defined by the complex vector z, as 

z = ±. (4.10) 

Let us remark that z = — J2k=i x k £ R. In the variable z, the potential takes the 
form 

1 N - x 1 
G^ N (z) = G^ N (x(Nz)) = -J2 A/mvN 2 + —\\x(Nz)\\i (4.11) 

k=0 



Moreover, by (1471) and (ROOD 



||x(JVz)||l = ||iVz||i^ = -i|M. (4.12) 

In the new coordinates the minima are now given by 

7± = ±(1,0,...,0). (4.13) 

In addition, z(B^) = z(B p ^(lJ)) = B p (I_) where the last ball is in the new 
coordinates. 

Lemma I47T1 will allow us to prove the following important estimates. For 5 > 0, we 
set 

C 8 = { z G R N : \z k \ < 5 ; fe ' JV , < k < N - 1 \ , (4.14) 



where X k:N are the eigenvalues of the Hessian at O as given in (12.141) and r k:N are 
constants that will be specified below. Using (I3.15D . we have, for 3 < k < N/2, 

r 2- 



X k , N >k z \l-—J^-l. (4.15) 
Thus (Xk,N) verifies X k)N > ak 2 , for 1 < k < N/2, with some a, independent of N. 
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The sequence (r k)N ) is constructed as follows. Choose an increasing sequence, 

{Pk)k>i> and set 

=1 ,<Kl»l (4 ' 16) 

rk,N = r N -k,N = Pk, 1 < k < Iy\ ■ 



Let, for p > 1, 

k'P 

\k>l 

Note that if K po is finite then, for all pi > p , K pi is finite. With this notation we 
have the following key estimate. 

Lemma 4.2. For all p > 2, there exist finite constants B v , such that, for z e C$, 

\\x{Nz)\\ p p <8»NB p (4.18) 

if K q is finite, with ± + ± = 1. 

Proq/: The Hausdorff-Young inequality (Lemma [4. ID gives us: 

\\x(Nz)\\ p <C g \\Nz\\ q ^. (4.19) 

Since 2; e C5, we get 

JV-l q 

\\Nz\\l T <8m^Y,T^- C4.20) 

fc=0 ^fe 

Then 

E^v 1 , 9 r r k,N 1 p fc 1 jyq = h p-n 

x g/2 X 9/2"^ Z x g/2 - ,,/2 a q/2 Lg - . g/2 "T a g/2 A 9 ^ l^-U 

fc=0 A fc A fc=l A fc A fc=l A 

which is finite if K q is finite. Therefore, 

\\x{Nz)\\l<5 p N {q - l)v «C p q D p q , (4.22) 
which gives us the result since (g — 1)| = 1. □ 
We have all what we need to estimate the capacity. 

4.1.2. Capacity Estimates. Let us now prove our main theorem. 

Proposition 4.3. There exists a constant A, such that, for all e < e and for all N, 

ca 1 p(B^B N )) . n-2 1 , , 

p ; +' ->> =eV2^s =(l + R(e,N)), (4.23) 

NN/2 -' v/|det(VF 7 ^(0))| 



where \R(e,N)\ < Ay/e\\ne\ 3 . 

The proof will be decomposed into two lemmata, one for the upper bound and 
the other for the lower bound. The proofs are quite different but follow the same 
idea. We have to estimate some integrals. We isolate a neighborhood around 
the point O of interest. We get an approximation of the potential on this neigh- 
borhood, we bound the remainder and we estimate the integral on the suitable 
neighborhood. 

In what follows, constants independent of N are denoted A*. 
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Upper bound. The first lemma we prove is the upper bound for Proposition 14.31 
Lemma 4.4. There exists a constant A such that for all e and for all N, 

-N-2 1 



j\W2-i - VZ 



TVS 



V|det(VF 7i7V (0))| 



(1 + A e\lne\ 2 ) . (4.24) 



Proof. This lemma is proved in [5] in the finite dimension setting. We use the same 
strategy, but here we take care to control the integrals appearing uniformly in the 
dimension. 

We will denote the quadratic approximation of C7 7i at by F , i.e. 



N-1 



N-1 



(4.25) 



fc=0 fc=l 

On Cs, we can control the non-quadratic part through Lemma [4721 



Lemma 4.5. There exists a constant A\ and 5 , such that for all N, 5 < 5 and all 

z e C s , 



G 7 ,n{z) - F (z 
Proof. Using (I4.11D . we see that 

G~f,N(z) - F (z) 



1 

iiV 1 



x(Nz)\\j. 



(4.26) 



(4.27) 



We choose a sequence {pk)k>i such that K±/ 3 is finite. 
Thus, it follows from Lemma [4721 with A 1 = \B±, that 



G 



7,iV 



N-1 



k,N\Zk\ 



k=0 



< A±5 4 



as desired. 



(4.28) 



□ 



We obtain the upper bound of Lemma 14.41 by choosing a test function h + . We 
change coordinates from x to z as explained in (14.10D . A simple calculation shows 
that 



\\Vh(x) 



N-'\\Vh(z)\\l 



(4.29) 



where h(z) = h(x(z)) under our coordinate change. 

For 5 sufficiently small, we can ensure that, for z ^ Cs with |^ | < 5, 



G,, N {z) > F (z) = , 2 



N-1 



fc,JV 



z k \ 2 > 



d 2 



+ 25 2 >5 2 . 



k=l 



Therefore, the strip 



Ss = {x\ x = x(Nz), \zq\ < 5} 



(4.30) 



(4.31) 



separates M N into two disjoint sets, one containing I and the other one containing 

I + , and for x e S d \ Cs, G^ iN (x) > 5 2 . 
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h + (z) = < 



The complement of S$ consists of two connected components T + , T_ which con- 
tain I + and J_, respectively. We define 

(1 for^er_ 

for zeT + , A . 

ft s r ' (4 - 32) 

f(z ) tor zeC 5 

arbitrary onSs\C$ but || Vh + \\ 2 < f. 

where / satisfies f(8) = and f(—S) = 1 and will be specified later. 

Taking into account the change of coordinates, the Dirichlet form (I2.9D evaluated 
on h + provides the upper bound 



e -G^M*)/z\\Vh+(z)\\ldz 



2((B^UB«) C ) 



e / e - G "<- N(z)/£ (f'(z Q )) dz + e5- 2 c 



(4.33) 

e -G j!N (z)/e dz 



The first term will give the dominant contribution. Let us focus on it first. We 
replace G 7i jv by F , using the bound (14.261) , and for suitably chosen 5, we obtain 

6 a> 



(f(z )Ydz < (1 + 2A, 



-F (z)/e ( f i 



(f (z )) dz 



^l + 2Ai^ J e _ A E ^=i lAfc - Jv|zfc|2 ^i...c/2jv_i 



x J (f'(z )] fe z o/ 2£ dz . 



(4.34) 



Here we have used that we can write C$ in the form [—8, 5] x D s . As we want 
to calculate an infimum, we choose a function / which minimizes the integral 

j 5 _ s (f'(z )) 2 e z o/ 2e dz . A simple computation leads to the choice 



Therefore 



-GLjv(*)/e / fi 



(f(z )Ydz< 



_ f 5 e- t2 l 2 Ht 
f{Z ° ] = J 7 s e-^dt 



L e~^^fe=o l A >=.ivllzfcl dz 

2, j ^ JCg 



S-5 e ^ 4 dzo) 



(4.35) 



5 4 

l + 2A 1 —\. (4.36) 



Choosing 5 = \/Ke \ \ne\, a simple calculation shows that there exists A 2 such that 

J c e ~^^=o\ x k,N\\z k \ 2 dz 



e 2 z o/ £ dz 



< V2 



-N—2 



TVS 



V|det(VF 7>7V (0))| 



[l + A 2 e). (4.37) 



The second term in Q4.33D is bounded above in the following lemma. 

Lemma 4.6. For 5 = \/Ke\ ln(e) | and p k = Ak a , with < a < 1/4, there exists 
A 3 < oo, such that for all N and < e < 1, 



A^^hTe 



N-2 



Ss\Cs ~ " #PUl 



3^/2+1 



(4.38) 
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Proof. Clearly, by (HTTTTh 



Thus 



Now, 



k,N\ z k\ 



k=l 



' Ss\Cs 
rS 

< I dz 



dzi . . . dz^-ie 2e 



1 EfeLo 1 ^k,N\ z kV 



5 J 3^=1 "-\ z k\>Srk,N I ' \J '^k,N 

6 N-l 



< 



[ e+^dzoT, [ 
J-s T~1 J\z, 



IT 



x I I f e- x - N ^ 2/2£ d Zi 

l<i^k<N-l J R 

N-l 




N-l 



\ i=l k=l 



N-l 

E r fc 

k=l 



L 2 



(4.39) 



(4.40) 



1 S 2 ri N /2e 



TJ N-l 
2_^t 1 k,N c ' 2—< 1 N-k,N c 



k=l 



k=l 



(4.41) 



We choose p k = 4k a with 0<a<l/4to ensure that K4/3 is finite. With our choice 
for 5, the sum in (14.41 D is then given by 



1 OO -.CO 



e 6Kn2a < Ce 2K , 



(4.42) 



n=l 



n=l 



since the sum over n is clearly convergent. Putting all the parts together, we get 
that 



-N-2 



ss\cs Vldet (V 2 F 7iJV (0))| 



and Lemma [4761 is proven 



(4.43) 

□ 



Finally, using (147361) , (147331) , (147371) , and (147381) , we obtain the upper bound 



N-2 



~ 7|det(VF 7iJV (0))| 



1 + A 2£ ) (1 + 2A 1 e\ \ne\ 2 + A 1 ^' 2 ) {AAA) 



with the choice p k = 4k a , < a < 1/4 and 5 = lne|. Note that all constants 

are independent of N. Thus Lemma 14741 is proven. □ 
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Lower Bound The idea here (as already used in fl6]]) is to get a lower bound by 
restricting the state space to a narrow corridor from I to I + that contains the 
relevant paths and along which the potential is well controlled. We will prove the 
following lemma. 

Lemma 4.7. There exists a constant A4 < oo such that for all e and for all N, 

1 



cap (B*,B?) r 



-N-2 



Tic 



y/\det(VF 1>N (0))\ 



1 - A A \ e\ lnel 3 . 



Proof. Given a sequence, (p k ) k >i, with r k . N defined as in (14.16D , we set 

C 5 = \z e] - 1 + p, 1 - p[, \z k \ < 5 r kiN /y/\ ktN \ . 



(4.45) 



(4.46) 



The restriction \z \ < 1 — p is made to ensure that Cs is disjoint from B± since in 
the new coordinates (14.101) I± = ±(1, 0, . . . , 0). 

Clearly, if h* is the minimizer of the Dirichlet form, then 



cap (5^, Sj) =*(/»*) >* e< (/»*), 
where $g is the Dirichlet form for the process on Cs, 



(4.47) 



-G lN (x)/e 



\Vh{x)\\ 2 2 dx = N N/2 - 1 e 



-G JiN (z)/e 



\\Vh{z)\\ 2 2 dz. 



To get our lower bound we now use simply that 

N-l 

\\vkz)\\t = Yl 



k=0 



dh* 


2 

> 


dh* 


dz k 




dz 



so that 



NN/2-l 



> e 



z{C s ) 



dh* 



dz n 



dz = $c s ( h *) - 



hen 



(4.48) 
(4.49) 

(4.50) 



The remaining variational problem involves only functions depending on the single 
coordinate z , with the other coordinates, z± = (^i)i<i<jv-i, appearing only as 
parameters. The corresponding minimizer is readily found explicitly as 

f !~P e G 1)N (s,z^)/e^ s 
J zn 



h (zq,z±) 



J 1 1 p e G 7 ,jv(s,2j_)/£^ s 



-1+p' 



and hence the capacity is bounded from below by 



cap {B N ,B^) 

NN/2-1 



>* d (h-)=e 



c 1 

A 



1-P - 



e G^ N (z ,z x )/e dz \ 
1+P ' 



-1 



dz i . 



(4.51) 



(4.52) 



Next, we have to evaluate the integrals in the r.h.s. above. The next lemma pro- 
vides a suitable approximation of the potential on Cs- Note that since z is no 
longer small, we only expand in the coordinates z±. 

Lemma 4.8. Let r k:N be chosen as before with p k = 4k a , < a < 1/4. Then there 
exists a constant, A 5 , and 8 > 0, such that, for all N and 5 < 5 , on Cs, 



N-l 



k=l 



< A 5 6 3 , 



(4.53) 
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where 



JV-1 



Zfcl 2 . 



(4.54) 



k=l 



Proof. We analyze the non-quadratic part of the potential on Cg, using (14. IIP and 
(1431) 

1 



N-l 



JV-1 



| i JV ~ 1 1 JV ~ 1 



i=0 



2=0 

>JV-1 



JV-1 



20 



fc=l 



JV-1 



i=0 



(4.55) 



where m = ^ u/ fc z fc . Note that V^" 1 u< = and u = ^x (JV(0, zj) and that 



20 

for z e R^, m,: is real. Thus 



JV-1 



JV-1 



^|l + Ml | 4 = iV+^(6^ + 4^ + M 4 ) 



(4.56) 



i=0 



we get that 



i=0 



i=0 



<|(4«+H!). 



A simple computation shows that 



N 



(4.57) 



(4.58) 



k^0 



Thus as \zq I < 1, we see that 



^||o;(A^)|| 4 -4-6z 2 £ 



fc^O 



<-(M\x(N(O jZ± ) 



Using again Lemma |4.2[ we get 

\\x(N(0,z ± ))\\l < B 3 N5 3 



\x(N(0,z ± ] 



< B A NS\ 



Therefore, Lemma [4781 is proved, with A 5 = 4£? 3 + B^5q. 
We use Lemma [4781 to obtain the upper bound 



\x(N(0,z ± ))\\t). (4.59) 



(4.60) 



□ 



/ ' e 5 ^ Z0 ^ )/e dz < ex P U e J2 k , \ ktN \z k \ 2 + ^f)g(z ± ), (4.61) 
J-i+p v J 



where 



9(z±) = J i ex P 1 Q*o - ^4 - 3o/(«J.))) ^o- 
This integral is readily estimate via Laplace's method as 

g(z ± ) = (1 + 0(e)) = V2^e{l + O(e) + 0(5 2 )) . 

V 1 - 2 f(z±) 



(4.62) 



(4.63) 
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Inserting this estimate into (I4.52D , it remains to carry out the integrals over the 
vertical coordinates which yields 

N-l 



$ g (/T) > el exp -l^A^N 2 - — ) (l + 0(e) + 0(6 2 )) dz ± 

^J dx (~ KMzkA dz ± (1 + 0(e) + 0(6 2 ) + 0(6 3 /e)) 




(4.64) 

The integral is readily bounded by 



N-l / 1 N-l 

- 5J / dzi . . . / . . . / ctejv-i exp -— V" A fc) 

fe=1 ■/» JM>Sr h ,N/y/>^N ■'K \ ^ £ fc=1 



zJ 2 ttej 



-jv-i -n- at i / /2e 



i=l V fe=l 



when 5 = VKehie and O(e^) uniform in N. Putting all estimates together, we 
arrive at the assertion of Lemma 14. 71 

□ 

4.2. Uniform estimate of the mass of the equilibrium potential. We will prove 
the following proposition. 

Proposition 4.9. There exists a constant A 6 such that, for all e < e and all N, 
where \R(N,e)\ < A 6 ^/e\ lne| 3 . 

Proq/: The predominant contribution to the integral comes from the minimum J_, 
since around J + the harmonic function h* BN BN (x) vanishes. 

The proof will go in two steps. We define the tube in the ,2 -direction, 

C s =^z: \/ k >i\zk\ < Sr ktN /y/\ htN } , (4.67) 

and show that the mass of the complement of this tube is negligible. In a second 
step we show that within that tube, only the neighborhood of I gives a relevant 
and indeed the desired contribution. The reason for splitting our estimates up in 
this way is that we have to use different ways to control the non-quadratic terms. 



Lemma 4.10. Let r k:N be chosen as before and let 5 = \jKe\ lne|. Then there exists 
a finite numerical constant, A 7 , such that for all N, 

1 f — 

/ e -GiM*V'dx < A 7 V2^ N ^ 64£ e K . (4.68) 

N N ' 2 Jci ~ Vdet(VF 7iJV (0)) 
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The same estimate holds for the integral over the complement of the set 

D & = {x : |z - 1| < 5 V |2o + 1| < 8}. (4.69) 
Proof. Recall that z = — J^iLo 1 x %- Then we can write 

G>rA*) = ~\ z l + \A + 2^(^ Da; ) - g^NIa + \ z o + -^W x Wt ~ \ z t- (4.70) 
Notice first that by applying the Cauchy- Schwartz inequality, it follows that 



4 = N~ 



'N-l 

E 

i=0 



1 / 



< N 



-2 



'N-l 

£ 

i=0 



N-l 



< 



(4.71) 



i=0 



Moreover, N 1 \\x\\l — \\z\\\, so that expressed in the variables z, 



N-l 



1 1 1 N ~ x 



2 u 4 

Therefore, as in the estimate (14.40D , 



4 



(4.72) 



k=l 



-G ltN (z)/e^ z <- 



N-l 



-e-\z%/A-z, 



k=l 



-Afc,jv|2 fc | 2 /2£^ 



' n 

l<i^k<N-l 



-\,n\zi\ 2 / 2e d z . 



(4.73) 



< 



< 



H*S/4-*g/2) d(Zo< j-l 



iV-1 



iV-1 



v ri 2 ^E^>~ 52r ^ /2e 

\ i=l fc=l 



- e -l(**/4-«g/2) 



AT-1 



\ n ^1^. 

\ i=l 



Since clearly, 



[ e - e_1 (*b 4 /4-«g/2) dzb = 2v ^e 1 / 4e (l + 0(e)), 



(4.74) 



this proves the first assertion of the lemma. 

Quite clearly, the same bounds will show that the contribution from the set where 
| z ± 1 1 > 8 are negligible, by just considering now the fact that the range of the 
integral over z is bounded away from the minima in the exponent. □ 

Finally, we want to compute the remaining part of the integral in (I4.66D , i.e. the 
integral over C s H {x : \z + 1| < 8}. Since the eigenvalues of the Hessian at J_, 
VkN> are comparable to the eigenvalues \ k N for k > 1 in the sense that there is a 
finite positive constant, c^, depending only on fi, such that X kjN < v k)N < c^A fcjv , 
and since u N = 2, this set is contained in C c g, where 



Cs(I-] 



z E R N : \z + l| < 



, W\ < 8 



rk,N 

y/Vk,N 



1< k <N- 1 



(4.75) 



UNIFORM ESTIMATES FOR METASTABLE TRANSITION TIMES IN A COUPLED BISTABLE SYSTEM 19 



It is easy to verify that on Cs(I-), there exists a constant, A 8 , s.t. 



N-l 2 

ri 



\z - < S 2 y < 5M 8 AT 2 2 . (4.76) 



and so, for 5 = \jKe\ lne|, Cg(I-) C z(B_). 

On Cs(I-) we have the following quadratic approximation. 

Lemma 4.11. For aZZ iV, 

AT-1 



G 7 .tf(*) + ^ - 2 ^kfel 2 = (4-77) 

k=0 

and there exists a constant A 9 and 5 such that, for 5 < 5 , on C«s(I_) 

\R(z)\<A 9 6 3 (4.78) 
where the constants (p k ) are chosen as before such that K 4 / 3 is finite. 

Proof. The proof goes in exactly the same way as in the previous cases and is left 
to the reader. □ 

With this estimate it is now obvious that 

h* BN BN (z)e- d ^ N{z)/£ dz = [ e- d ~" N{z)/£ dz (4.79) 

■N 



c l/4e 



v /detV 2 F 7iA r(/_ 



(1 + 0(5 3 /e)) 



Using that h* BN b n(z) vanishes on and hence on Cs(I+), this estimate together 
with Lemma [4. 101 proves the proposition. □ 

4.3. Proof of Theorem [37T1 

Proof. The proof of Theorem I3TT1 is now an obvious consequence of (I2.16D together 
with Propositions [473] and [4791 □ 
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